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We study quantum trajectories of collective atomic spin states of N effective two-level atoms driven 
with laser and cavity fields. We show that interesting "entangled-state cycles" arise probabilistically 
when the (Raman) transition rates between the two atomic levels are set equal. For odd (even) N, 
there are (N + 1) /2 (N/2) possible cycles. During each cycle the 7V-qubit state switches, with each 
cavity photon emission, between the states (\N/2, m) ± \N /2, —m))/\2, where \N/2, m) is a Dicke 
state in a rotated collective basis. The quantum number m (> 0), which distinguishes the particular 
cycle, is determined by the photon counting record and varies randomly from one trajectory to the 
next. For even N it is also possible, under the same conditions, to prepare probabilistically (but 
in steady state) the Dicke state |iV/2, 0), i.e., an iV-qubit state with N/2 excitations, which is of 
particular interest in the context of multipartite entanglement. 

PACS numbers: 42.50.Dv, 42.50.Lc, 42.50.Pq 



I. INTRODUCTION 

The understanding that entanglement serves as a phys- 
ical resource for various quantum communication or in- 
formation processing protocols has motivated both ex- 
perimental and theoretical studies in the preparation and 
quantification of entangled states [TJ [2J. Entanglement 
between two subsystems of a quantum system, i.e., bipar- 
tite entanglement, is well understood and can be directly 
quantified by a number of readily computable measures. 
Multipartite quantum systems, however, offer a richer va- 
riety of entangled states, the characterization of which is, 
naturally, more complicated and indeed still under active 
investigation (see, e.g., [1 H EJ EJ [S] ) • 

Of particular interest in this context are the A-qubit 
symmetric Dicke states [9]; in particular, the states \j, m) 
(with j — N/2) satisfying 

Jz\j,m) = m\j,m), 3 2 \j,m) = j(j + l)\j,m), (1) 

with m £ {—j, ■ ■ ■ where { J z , J 2 } are col- 

lective angular momentum operators and (m+N/2) gives 
the total number of excitations in the system. These 
states are, with regards to entanglement, particularly ro- 
bust against particle loss [7], and they are significant for 
the study and application of genuine multipartite entan- 
glement OH EE]. 

Photonic systems have proven fruitful for the ex- 
perimental preparation and detection of three-photon 
polarization-entangled W states (Dicke states with one 
"excitation" ) jT2J [T31 [21 [T5J [TB] and, very recently, four- 
photon polarization-entangled Dicke states with two ex- 
citations [T7j. Meanwhile, single collective atomic ex- 
citations - atomic W states - have been produced in 
many-atom ensembles [T5J [HJJ [2U1 123 H2] and with small 
collections of trapped ions [25] • The potentially long co- 
herence times of ground electronic states of atoms (two 
of which constitute the qubit) make atomic systems par- 
ticularly attractive for the storage and manipulation of 



entangled states, and a variety of schemes have been pro- 
posed for the preparation of more general Dicke states of 
trapped atoms or ions [24l 125] 126 ] l27 j 128 ] 129 ] [30] . 

Here we describe a system that adds to these proposed 
schemes, in that it may produce collective atomic states 
with multiple excitations, but also yields very interest- 
ing "entangled-state dynamics" and conditional quantum 
evolution, based upon photon detection in the output 
field of an optical cavity containing the atoms. The cav- 
ity QED setup and excitation scheme that we employ 
have been considered previously in the context of gener- 
ating, in steady state, two-qubit (i.e., two-atom) entan- 
gled mixed states of any allowed combination of purity 
and entanglement [31] , although here we consider a spe- 
cific limiting case of the operating conditions (not consid- 
ered in |31j ) that does not in fact admit a unique steady 
state. 

In examining this specific case, we gain inspiration 
from a recent investigation of conditional quantum dy- 
namics in a related cavity QED system involving two 
atoms in separate, cascaded cavities [52] , which also fo- 
cussed on a specific limiting case of the operating condi- 
tions for an atom-entanglement scheme proposed earlier 
[33] . In particular, we use the method of quantum tra- 
jectories based upon continuous monitoring of the cavity 
output field by photodetection to study the conditional 
quantum evolution of the collective state of N atoms in- 
side the (single) cavity. As in [32], we observe manifestly 
distinct behaviours from one trajectory to the next; the 
system either evolves to a (steady) Dicke state with N/2 
excitations without emitting any photons (for even N), or 
it executes a sustained cycle of "switches" between par- 
ticular, well-defined superpositions of Dicke states, with a 
concomitant continuous output stream of photons. Fur- 
thermore, for N > 2 more than one distinct switching 
cycle is possible, with the different cycles distinguished 
by the rate of photon emissions. 

We start in Section [IT] with a brief review of the cav- 



2 



ity QED configuration proposed in [51] , after which we 
formulate our quantum trajectory model, making use of 
a convenient change of basis, which is motivated by the 
particular form of the effective Hamiltonian and jump 



operators. In Section III we focus briefly on the particu- 
lar case in which the system evolves into a stable Dicke 
state with N/2 excitations. Individual trajectory results 
for N — 2, 3, 4 atoms, highlighting cyclic behavior of the 



state, are presented and discussed in Section IV We also 



make some more general observations concerning cycles 
for larger N , before concluding in Section W\ 
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II. PHYSICAL SYSTEM AND THEORETICAL 
MODEL 

A. Multiatom cavity QED system 

We consider a system of N four-level atoms trapped 
inside a high finesse optical cavity. The atomic excita- 
tion scheme and physical implementation are depicted in 
Fig. [I] A single quantized cavity mode couples to the 
|0) <-* |r) and |1) «-> \s) transitions in each atom with 
coupling strengths g r and g s , respectively. Auxiliary laser 
fields couple to the |0) \s) and |1) «-> \r) transitions 
with Rabi frequencies £ r and ( s , respectively. Together, 
the cavity and laser fields drive (resonant) Raman tran- 
sitions between the stable atomic ground states |0) and 
|1). Such an excitation scheme could be achieved, for ex- 
ample, by using orthogonal, linearly-polarized laser and 
cavity fields on an F = 1/2 <-* F' = 1/2 atomic transi- 
tion (as in 6 Li), or on an F = 1 <-> F' = transition (as 
in 87 Rb) with an applied magnetic field, as described in 
[33). 

Finally, we note that the spacings between the atoms 
are assumed to be sufficiently large that direct dipole- 
dipole interactions can be neglected. We also assume 
that each atom is held tightly in position by a confining 
trap, so that the atomic motion can be omitted from the 
analysis. 



B. Reduced iV-atom master equation 

With the assumption of large detunings A r and A s , the 
atomic excited states can be adiabatically eliminated and 
atomic spontaneous emission neglected. The coherent 
dynamics is then characterized by Raman transition rates 
and energy level shifts, which we subsequently assume 
to be much smaller than the cavity field decay rate, k. 
This enables us to adiabatically eliminate the cavity field 
to derive a master equation that describes the atomic 
dynamics alone. Under "resonant" conditions (i.e., light 
shifts of the atomic ground states are balanced, such that 
Raman resonance is maintained; see |31j for details), this 
master equation can be written in the form 



dp 



2XpX^ - X^Xp - pX^X, 
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FIG. 1: (a) Atomic excitation scheme and (b) physical imple- 
mentation of the system studied. A single cavity field mode 
(coupling strength g r ,s) and two laser fields (Rabi frequencies 
( r . s ) drive resonant Raman transitions between the atomic 
ground states |0) and |1). Note that an additional laser 
or magnetic field may be required in the excitation scheme 
to provide an independent relative energy shift between the 
ground states |0) and |1) (see [31]). Note also that \r) and |s) 
may be a single level, provided the two Raman channels re- 
main distinct from each other (which would require that the 
two ground states are nondegenerate in energy). For large 
detunings A r and A s , and for k, the cavity field decay rate, 
much larger than the Raman transition rates, the model may 
be reduced to that of an effective two-level atom with char- 
acteristic transition rates 701 and 710. 



where p is the atomic density matrix, 

X = VToi J+ + VTio J- (3) 
with 701 = |Cs.9s/2A s | 2 /«, 710 = \( r g r /2A r \ 2 /«, and 



N 



En) 



j = (j + )t. 



(4) 



i=i 



For 701 7^ 710, the dynamics of ^ produces unique 
steady states that are non-maximally-entangled (pure or 
mixed) collective spin states [311 [35l [36] . However, for 
the remainder of this paper we consider the case 

7oi = 7io , (5) 

for which ^ does not admit a unique steady state (math- 
ematically, the system of linear equations for the density 
matrix elements possesses a second zero eigenvalue when 
7oi = 7io) ■ Scaling time by (4701) -1 , the master equa- 
tion we thus focus on takes the form 

dp 
~dt 



^"JxP'J X 



(6) 
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where J x — (J+ + J_)/2 is the ^-component of the to- 
tal pseudo-spin J of all TV qubits. Note that since J 2 
commutes with J x , the dynamics described by ([6| con- 
serves the magnitude of the total spin angular momentum 
(which we will take to be the maximal value in this work) 

C. Quantum trajectories in the collective-spin 
representation 

1. Change of basis 

We are interested in the time-dependent behavior of 
the atomic state subject to dynamics as described by ^ 
and, as depicted in Fig. [ljb), in photon counting mea- 
surements on the cavity output field. The method of 
quantum trajectories (see, e.g., [38j [39l HOI SI] ) , summa- 
rized briefly in the Appendix, lends itself nicely to this 
task. 

In this approach, and with photon counting measure- 
ments in mind, we rewrite ([6| in the form 

J = -iHctip + ipHts + 2Cptf , (7) 

and identify a non-Hermitian effective Hamiltonian, Ti. c s , 
and collapse operator, C, as 

H e s = —iJx , C — J x . (8) 

In the trajectory model, free evolution of the atomic 
wave function, subject to H c b, is interspersed, at ran- 
dom times, with collapses (or quantum jumps) produced 
by the action of C on the wave function. These collapses 
can be identified with photon detections in the cavity 
output field and the resulting wave function trajectory 
describes the state of the atomic spin system conditioned 
upon the sequence of photon detections. 



Given the form of both 7i e g and C, a particularly con- 
venient basis for our study will be the basis of eigenstates 
of the ^-component of the total spin J; more precisely, 
simultaneous eigenstates of J 2 and J x . In the more fa- 
miliar representation of simultaneous eigenstates of J 2 
and J z , which we denote by |j, m) z , the allowed values of 
the J 2 -component are m — —j, —j + 1, . . . ,j, where, for 
a system of N effective spin- 1/2 particles, as we consider 
here, j may take one of the values {0, 1/2, 1, . . . , N/2}. 
So, for example, for N = 2 atoms and m — 0, 

|0,0), = (|10), - |01)J , (9a) 
|l,0) z = -L(|10) Z +|01>J , (9b) 

while for m — ±1 (both spins "up" or "down") 

|1,-1) 2 = |00) Z , |1,1) 2 = |H>* , (9c) 

where \ij), = 0, 1) denotes the state with atom 1 in 
state \i) and atom 2 in state \j) . 

To obtain the eigenstates of J Xl we perform the ap- 
propriate set of Euler rotations on the coordinate sys- 
tem. In particular, the J^-eigenstates are related to the 
Jz -eigenstates by a unitary rotation operator with the 
property £> _1 (/3) = D(-/3): 

\j,n) z = J2 D H(-l)\j,m) x . (10) 

m=-j 

Here, D%> n (0) are the matrix elements of the rotation 
operator given by (sometimes referred to as Wigner's for- 
mula) gg 



j+m 

(/?) = £(- 

k=0 



■1) 



k—n-\-m 



\J U + n)\ (j - n)\ (j + m)\ (j - m)\ 



k\(j 



k)l(j 



k)\ (fe- 



ces ■ 



2j—2k-\-n—m 



2k— i i-\-rn 



(ii) 



Note that the sum in (11) is only over terms with non- 



negative arguments of the factorials. 

Finally, we will confine our study to evolution from the 
initial state 



|^(0)) = |7V/2,-iV/2) z = |0 



(12) 



i.e., j — N/2 (the maximal value) and all spins "down". 
Since j is conserved by the evolution, only N + 1 lin- 
early independent basis states are required to describe 
the state of the TV-atom system. 



2. Jumps versus free evolution: general behavior 



In the basis of Jaj-eigenstates it is very straightforward 
and instructive to calculate the general behavior of the 
spin wave function either undergoing free evolution sub- 
ject to HcS, or being "collapsed" under the action of the 
jump operator C — J x . In particular, consider first free 
evolution from the initial state 1-0^(0)) = \N/2,—N/2) z 
up to some time t, without any collapses having occurred. 
For even N, the (unnormalized) state IV'ivW) can b e 
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written as 

N/2 
m=-N/2 

(13) 

Noting that D^% 2 (-f ) = D^l N /2 K ), * is use- 
ful to define 

| X ± (m)) = ( |JV/2, m) a ± \N/2, -m) x ) , (14) 
and rewrite (l3| as 

|^(t)) = ^ /2 (-|) liV/2,0), 

TV/2 

+ E C ™^" l2t !^( TO ))' ( 15 ) 
m— 1 

where c m = V^D 1 ^^^ (— §)• This result shows explic- 
itly that the to = component of the spin wave function 
in the J^-eigenstate basis is unaffected by the free evo- 
lution, while the amplitudes of the to ^ components 
decay exponentially with time. 

Conversely, the m = component is completely 
removed by the action of the jump operator, i.e., 
J x \N/2,0) x = 0, while 

Jx\x%(m)) = m\x^(m)) . (16) 

For odd N, there is no to = component and the 
(unnormalized) wave function takes the form 

N/2 

hM*)>= E c m e- m2t \ X +(m)), (17) 

m=l/2 

after a time t of free evolution (from the initial state 
|V>w(Q)) = \N/2,-N/2) z ). Hence, all components of the 
wave function have amplitudes that decay exponentially. 

Now, in the trajectory formalism, the square of the 
norm of the wave function, ( , 0jv(t)|i/'Af(i)) < 1 ; can be 
identified as the probability for a collapse (i.e., a photon 
detection) not to have happened [43 . From this identi- 
fication and our observations above, it follows that the 
system under investigation can exhibit very distinct be- 
haviors, as we detail in the following two sections. 



that no collapse ever occurs and, hence, that the system 
evolves to a stable steady state given by |7V/2, 0)^ . This 
is a Dicke state, in the basis of J^-eigenstates, with N/2 
excitations, which, as mentioned in the introduction, is of 
particular interest in the context of multipartite quantum 
entanglement. In fact, for our system a necessary and 
sufficient condition for genuine TV-partite entanglement 

is [Tnum 

{Jl) < f , (19) 

which is "optimally" satisfied for the state |iV/2, 0)^,, for 
which (J x ) = 0. Note that with a standard unitary 
rotation, implemented for example using coherent laser 
Raman pulses, this state can be transformed into the 
corresponding state |iV/2,0) 2 in the (conventional) J z - 
eigenstate basis. 

The probability P;v(m = 0) is plotted in Fig. [2] as a 
function of N and maintains a significant value even for 
quite large N. In fact, using Stirling's approximation 
(n! ~ 7i™e n v / 2'7rn for large n), one can show that for 
iV>l, 

*M™ = °)^V^' (20) 

which describes the slow decrease with TV, as illustrated 
in Fig. [| 

It must be emphasized that in the case that the state 
\N/2,0) X is prepared, no photons are emitted from the 
cavity. With a probability 1 — Pjv("^ = 0) there is a 
collapse, or photon detection. This in fact initiates a 
sustained sequence of collapses (and photon detections), 
since the amplitudes of all of the remaining (to ^ 0) com- 
ponents of the wave function decay exponentially under 
free evolution, i.e., the norm of the wave function is now 
guaranteed to decay towards zero, meaning that a col- 
lapse is also guaranteed to occur. Hence, the system 
exhibits a distinctly bimodal behavior: either no pho- 
tons are detected, "signalling" preparation of the state 
| N/2, 0} x , or a continuous stream of photons is detected. 
While this stream of photons indicates failure in the 
preparation of a steady state, it also signals the initia- 
tion of some remarkable dynamics of the system, which 
we label entangled state cycles and examine in the next 
section. 



IV. ENTANGLED-STATE CYCLES 



III. 



STEADY STATE PREPARATION OF THE 
DICKE STATE liV/2,0)^ 



From (151, it follows that for an even number N of 



atoms there is a finite probability, given by 

TV! 

~ 2 N [(JV/2)!] 



Pjv(m = 0) = 



n (N/2) 



(18) 



As described above, a first photon detection initiates a 
persistent sequence of photon detections and the relevant 
atomic spin states involved in the dynamics from then 
on are the states |xaK to ))- Note that for odd N, a first 
photon detection and subsequent persistent sequence is 
guaranteed to occur. 

The rate of photon detections is proportional to 
the expectation value of the operator C^C = J x , 
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FIG. 2: Probability Pjv(m = 0) for preparation of the Dicke 
state \N/2,Q) versus N/2 (N even). The crosses are calcu- 
lated from Eq. ^lsj. The solid line is the curve [2/ (ttN)] 1/2 , 
which approximates Pjv(m = 0) for N 3> 1. 
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FIG. 3: Example trajectory for the initial state [1,— 1) = 
|00) z . Time is measured in units of (47io) _1 . No collapses 
occur and the system is eventually projected into the 11,0}^. 
eigenstate. 



and hence to the value to 2 . The photon detection 
configuration we consider can therefore be interpreted 
as providing a measurement of m , so we can ex- 
pect each trajectory to "converge" towards a particu- 
lar pair of states, |XaK to )) ( or possibly to the state 
|iV/2,0) for even N), conditioned upon the measure- 
ment record. For example, from ( 16 1 we see that collapses 
increase the amplitude (xjv( to + 1) I V'iv(i)) relative to 
{x%( m ) \ V'iv(i)), while from (15), we see that free evolu- 
tion decreases (Xjv( r 
to (x%(m)\^ N (t)) 



1) ipi\[(t)} exponentially relative 
Hence, within a given time interval, 
the more photon detections that occur, the more proba- 
ble it is for the state to converge towards states \x%( m )) 
with larger to (for which, of course, a higher photon count 
rate is also expected). In contrast, the longer \i[>N(t)} 
evolves without a collapse the more likely it is for the 
state to converge towards states |x]y( TO )) with lower val- 
ues of TO. 

Which state or pair of states (\N/2,0) x or |XaK to ))) is 
"selected" in a particular trajectory is random. In the 
case of |x A r( m ))' the sequence of photon detections cor- 
responds to a sequence of jumps, described by (16 1, be- 
tween the states |xjv( m )) an d \xjf( m ))- These states are 
also very interesting entangled states, and we label the 
system's behavior in this instance as an entangled- state 
cycle. We now look at specific examples of entangled 
state cycles for N — 2,3, 4. 



A. Entangled-state cycles for JV= 2 

Consider free evolution of the initial state 

|Va(0)) = |1,-1>„ = ^ (11,0), + |x 2 + (l)» (21) 



from t = up to time t, without any collapses having 
occurred. The normalized wave function is given by 



\Mt)} 



-iHesst 



IV> 2 (0)} _ 11,0), 



,-t\^+ 



xi(i)) 



||e-*«-*|^(0))|| 



(22) 



This illustrates very clearly that as time progresses and 
no collapses occur (i.e., no photons are detected), the 
more likely it is for the system to be in the state |l,0) x . 
Such a no-collapse trajectory (an example of which is 
shown in Fig. [3]) occurs with probability 



P 2 (m = 0) = 1(^(0)11,0^1 



2 ' 



(23) 



i.e., in 50% of trials the system is prepared in the two- 
qubit maximally-entangled state 



|1>0) 3 



1 



V2 



11,-D 



IM>J 



l 



^/2 



(|00L-|H)J 



(24) 



When a collapse does occur and a photon is detected, 
the system state is projected onto 



|xa(l)> = 



c\Mt)) = JL 

\C\Mt))\\ V2 



(IM^-IWU 



(25) 



With the m = component removed, subsequent time 
evolution is guaranteed to lead to further collapses and 
a persistent cycle is established between the two states 
|X2(1))- Such a cycle is illustrated in Fig. [Ij where the 
abrupt changes in the sign of the amplitude (V^ll, - l)a; 
indicate jumps between the states |x 2 (l))- 
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FIG. 4: Example trajectory for the initial state |1, — 1) = 
|00} z . The system settles into an entangled-state cycle after 
the first collapse. 



It is interesting to consider the compositions of the 
states |x2 (1)) m the original J 2 -eigenstate basis; in par- 
ticular, 

|X2(1)> = |1>°>*=4|(I 1>* + |1 >J > ( 26a ) 

and 

|x 2 + (i)) = ^ (|i, i), + |i, -l) J - ^ (100), + |11) J . 

(26b) 

So, when a cycle is established, the states involved are 
also maximally-entangled states of the two qubits. 



B. Entangled-state cycles for JV=3 

Considering the initial state 



|V*(0)) = |3/2,-3/2) 
ft 



1 



(27) 



we see that there is no m = component and hence an 
entangled state cycle always results. In particular, a first 
collapse is guaranteed, as exemplified by the wave func- 
tion after time t of free evolution (ignoring normalization 
for simplicity, i.e., Eq. (17 1 for N = 3), 



l-v+a 



With probability 

l(X3 + (^)l^3(0))| 2 
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(29) 
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FIG. 5: Example trajectory for the initial state 
3/2, — 3/2> 2 = |000} z . The m = ±3/2 components de- 
cay to a negligible magnitude by the time of the first collapse 
and the system evolves into the cycle between the states 
|x 3 ± (l/2)). 



the system evolves to an entangled-state cycle between 
the states |xf (1/2)). A sample trajectory of this cycle is 
shown in Fig. [5] In this particular example, by the time 
of the first collapse, or photon detection, the components 
of |3/2, ±3/2)^ have become negligibly small. 

A trajectory in which the three-qubit state evolves to a 
cycle between the states (3/2)) (occurring with prob- 
ability 1 /4) is shown in Fig. [6 In this case, collapses are 
clearly much more frequent than in the previous exam- 
ple, since, mathematically, the norm of the wave function 
decays more rapidly between each collapse. In fact, the 
ratio of the frequencies of collapses for the two cycles is 
given by 





etc 


xUD) 


<xf(l) 


etc 


xf(D) 



(30) 



so which cycle the system evolves into for any individ- 
ual realization is clearly distinguishable from the photon 
counting record. 

It is interesting to consider the states involved in the 
two cycles. In the J^-eigenstate and J z -eigenstate bases, 
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FIG. 6: Example trajectory for the initial state 
3/2,— 3/2) = 1 000) , where sufficient jumps occur over a 
time interval to promote the three-qubit state into a cycle 
involving jx^(3/2)). 



respectively, we have 



1 



|x^(l)> = ^[(|no) x + |ioi) :c + |oii) x ) 

+ (1001), + 1010)^ + 1100) J] (31a) 



2V5 



(|011),+ |101), + |110),) 



73 



1000), 



(31b) 



k 3 -(|)> = -^[(|no) :c + |ioi) x + |oiiy 



(|001) X + |010) S +|100)J] (31c) 



2^3 



(|10Q) +|01G), + |001),) 



73 



(31d) 



and 



\xtm = -^(iii^.+iooo) j 



1 

72 
1 

2 

1 



(32a) 

(|000), + |011), + |101), + |110),) , 



|xs (i)> = -7| (ini)x -loooy 



(32b) 
(32c) 



= ~(|111), + |100), + |010), + |001),) . 

(32d) 

In the J x -eigenstate basis, the states \xf (3/2)) are sim- 
ply GHZ states, while \xf (1/2)) are each a linear com- 
bination of two W states, |3/2,l/2) and |3/2,-l/2) . 



Such combinations are sometimes referred to as G states 
[45 , which fall under the GHZ class of entangled states 
[3J, but are clearly less entangled than a single W 
state. Alternatively, in the J z -eigenstate basis, the states 
|x 3 (1/2)) and |x3 (3/2)) are all superpositions of aprod- 
uct state and a W state (T2j [14] . Note again that uni- 
tary rotations of states between bases could be imple- 
mented using coherent Raman transitions between the 
atomic ground states |0) and |1). 



C. Entangled-state cycles for JV=4 

In the case of four qubits with initial state |2, —2). 
1 0000) ,, pi reads 



\Mt)) = 



i 

71 



|X4 + (1)) 



12,0), 



1 

71 



|x 4 + (2)) 



(33) 



so again two distinct entangled-state cycles are possi- 
ble, corresponding to m — 1 or m — 2 and occur- 
ring with probabilities 1/2 or 1/8, respectively. In the 
J x -eigenstate basis, the cycles occur again either be- 
tween two GHZ states, |x4 t (2)), or between two G states, 

|xJ(i)>- 

However, when represented in the J z -eigenstate basis, 
the cycles now occur between quite different forms of en- 
tangled states. In particular, the m = 1 cycle is between 
a GHZ state and a G state, i.e., 



1x1(1)) 



^ (loooo) , 



78 



[(10001), 
11110), 4 



-lllll)J , 

+ |0010), 
11101), 4 



+ |0100) . 
11011), 4 



(34a) 

+ |100Q)J 

|0H1)J] , 
(34b) 



while the m = 2 cycle is between a G state ( XT (2))) and 
a state that is an equal superposition of a GHZ state and 
a Dicke state with two excitations [10] [T7], i.e., 

\xt{2)) = -5= (|1100) 2 + |1010), + |1001), + |0101), 



|0011), + |0110), + |0000) z 



|XI(2)) 



78 



(10001), 



|1110) 2 



+ |0010), 
11101).+ 



+ |0100), 
11011). 4 



limy , 

(35a) 
+ |1000), 

|om)j . 

(35b) 



(Note that the state |2, 0) x is of the same form as ( |35a I , 
but with different weightings for the GHZ and Dicke 
states.) Examples of entangled-state cycles for N = 4 
are shown in Figs. [7] and [8] 
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FIG. 7: Example trajectory for N = 4 with initial state 
2, — 2} = 1 0000) ^ . In this case, the system settles into the 
cycle between the states \xtW)- 



<m o 



10 



10 



10 



10 



<m" 



ii'iiii 



4 6 
time 



10 



FIG. 8: Example trajectory for N = 4 with initial state 
2,— 2) z — 0000) z . In this case, the system settles into the 
cycle between the states |xf (2)). 
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FIG. 9: Probability P(XN( m )) f° r the system to evolve into 
the entangled-state cycle between |Xjv( m ))- The crosses are 
calculated from Eq. ( 36 1 The solid lines are the Gaussian ap- 



proximation given in Eq. ( 37 1 



D. Larger N 

For larger N the number of possible cycles obviously 
increases, as does the variety of entangled states involved. 
From the decomposition of the initial state |0® ) z in 
terms of the states \N/2, 0) x and |xjv( m )); one deduces 
that the probability, P(x%{m)), for the system to be 
projected into the cycle involving the states \x%(m)^ is 
given by 

P(x N (m)) = 2 



D 



{N/2) 
m,~N/2 



TV 

2 
N\ 



2 N - 1 (N/2 + m)\{N/2 - to)! 



(36) 



In the large-TV limit, and for to <C N, this is approxi- 
mated well by the simple Gaussian form (derived using 
Stirling's approximation) 



P(X N (m)) 



exp 



to 2 \ 
N/2) 



(37) 



as demonstrated in Fig. [9] for N = 20 and N = 40. This 
result makes it very clear that cycles with m > \fN be- 
come increasingly unlikely for large N . Favored cycles 
are those with small to, corresponding, in fact, to states 
satisfying the condition for genuine iV-partite entangle- 
ment, 



{Jl) < 



N 



(38) 



i.e., states with (J%) = to 2 < AT/4. 
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V. CONCLUSION 

In this work we have described a means of realizing, 
through suitable laser and cavity QED interactions, a 
collective atomic spin system with dynamics described 
by the simple master equation 



d P O 7 7 72 
— = 2J x pJ x - J x p 



Despite its apparent simplicity, this model reveals a re- 
markable diversity of behavior when viewed in a quantum 
trajectory picture, given an initial state | AT/2, —N/2) z — 
\0® N ) Z and assuming photon counting measurements on 
the cavity output field. The different behaviors can be 
readily understood by employing a basis of J^-eigenstates 
and interpreting the configuration in terms of a projective 
measurement of the operator J x . On top of the inter- 
esting dynamical behavior, the actual states produced, 
either in steady state (m = 0) or participating in an 
entangled-state cycle (m > 0), typically possess genuine 
multipartite entanglement. 

Our work is actually quite closely related, but com- 
plementary, to that of Stockton et al. |28j . who con- 
sider continuous quantum nondemolition measurement of 
the atomic spin operator J z , via a dispersive atom-cavity 
interaction and homodyne detection of the cavity field, 
starting from an initial atomic spin state polarized along 
the x-axis. There, Dicke states in the J z -eigenstate basis 
are prepared, cither nondctcrministically, with a prob- 
ability determined by the initial state decomposition in 
this basis, or deterministically, via the use of state-based 
feedback control. Such a feedback approach could per- 
haps also be applied usefully to the present configuration. 
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APPENDIX A: MONTE CARLO 
WAVE-FUNCTION ROUTINE 

We will not discuss quantum trajectories in any pro- 
found way, but we simply outline how a realization of 



the reduced density operator for the system, called a tra- 
jectory, can be built operationally using a Monte Carlo 
algorithm. In essence, the trajectory formalism reduces 
the problem of mixed-state evolution to that of a pure 
state, as we describe below. References [3H U01 E] 
provide an excellent account of the formalism. 

For our quantum trajectory approach, the output of 
the optical cavity is considered to be monitored continu- 
ously by photon counting, as shown in Fig. |TJb) . In this 
scenario, each photodetection sequence corresponds to a 
trajectory of the A^-qubit collective atomic state. 

Consider a normalized pure state \ip(t)) at time t. The 
state is evolved "freely" over a small (ideally infinitesi- 
mal) time increment dt according to 



\ip(t + dt)) = 



(i-iH cS dt/h)\i>(t)) 
\{i-iu cS dt/h)\^{t))\y 



(Al) 



where 7Y c ff is a non-Hcrmitian effective Hamiltonian and 
the state is renormalized after the non-unitary evolution 
|46j . Such free evolution is interspersed, at random times, 
with "quantum jumps" or "collapses" , the probability of 
which must be calculated at each timestep. In particular, 
the probability that a collapse of the state occurs over the 
(infinitesimal) time step dt is given by 



v{t) = m)\&c\ii>{t))dt, 



(A2) 



where C is the collapse operator. To determine whether 
a collapse occurs or not, a uniformly distributed random 
number u(t) is drawn from the unit interval [0, 1] and 
compared to V(t). If V(t) > u(t) then the state becomes 



\ip(t + dt)) = 



c\m) 
\\c\m)\\ 



(A3) 



otherwise, for V(t) < u(t), the state is not collapsed 
and it continues to evolve, according to (All, over an- 



other subsequent interval of dt. This process is iterated 
over sufficiently many time steps to build up a trajec- 
tory. Note that for our system the action of C represents 
the emission and subsequent detection of a photon in the 
cavity output field. 
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